Low-temperature crossover in the momentum distribution of cold atomic gases in one dimension 
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The momentum distribution function for the two-component ID gases of bosons and fermions is 
studied in the limit of strong interatomic repulsion. A pronounced reconstruction of the distribution 
is found at a temperature much smaller than the Fermi temperature. This new temperature scale, 
which equals the Fermi temperature divided by the dimensionless coupling strength, is a feature of 
the two-component model and does not exist in the one-component case. We estimate the parameters 
relevant for the experimental observation of the crossover effect. 



I. INTRODUCTION 

During the last few years the great advances in the art 
of cooHng and trapping atomic gases have made it possi- 
ble to study experimentally the behavior of fermion and 
boson gases in lower dimensions Of particular interest 
are the quasi-one-dimensional systems which can be real- 
ized as one-dimensional (ID) tubes in a two-dimensional 
optical lattice j2, IS Q- The ground states of such 
systems are characterized by the ratio between the inter- 
action energy and the kinetic energy. At high densities, 
the kinetic energy dominates and the system can be de- 
scribed by mean-field theory. At low densities, one enters 
a regime of strong coupling where the interaction energy 
dominates the kinetic energy. The difference between this 
strongly interacting regime and the regime of weak inter- 
actions is more than just quantitative. One should expect 
that strong interactions result in phenomena absent in 
the weakly interacting case. The question of what these 
phenomena are and their experimental consequences is 
essential for our understanding of low-dimensional sys- 
tems and, more generally, the physics of strongly inter- 
acting matter. 

It is evident that the presence of internal degrees of 
freedom enriches the phase diagram of a physical system 
and may even give rise to completely new phenomena. 
One-dimensional systems in general and one-dimensional 
quantum gases in particular, are no exceptions to this 
rule. In this paper we discuss a phenomenon in quantum 
gases which is specific to the strongly interacting regime 
and the presence of internal degrees of freedom. We show 
that, due to the emergence of a soft propagating mode 
in the limit of strong repulsion, one-dimensional gases 
(both bosonic and fermionic) undergo a dramatic recon- 
struction of the one-particle momentum distribution as 
the temperature is varied on a very small scale. This tem- 
perature scale is proportional to the propagation velocity 
of the soft mode and vanishes in the limit of infinitely 
strong interparticle repulsion. 

The paper is organized as follows. In Section ^ we 
give a qualitative comparison of one- and two-component 
quantum gases in ID and introduce the relevant energy 
scales. We show that in the two-component case there 
exists a new low-temperature regime where the behavior 
of the quantum gas should be qualitatively different from 



that at zero temperature. Sections ITTTI and ITvl are techni- 
cal. In Section ITTTI we present an exactly solvable model, 
which grasps the physics discussed in Section |n| We give 
an exact representation for the one-particle density ma- 
trix in this model. In Section}^ we analyze short- and 
long-distance asymptotics of the density matrix and dis- 
cuss the numerical recipe for evaluating it at intermediate 
distances. Section Fvl contains the main result of the pa- 
per: The effect of the low temperature reconstruction of 
the momentum distribution function is demonstrated for 
both boson and fermion gases. In Section IVTI we discuss 
the possibility of observing this reconstruction effect ex- 
perimentally and estimate the magnitude of the coupHng 
strength for such an experiment in terms of the parame- 
ters characterizing highly elongated harmonic traps. 



II. ONE-COMPONENT VERSUS 
TWO-COMPONENT GASES IN ONE 
DIMENSION 

One-dimensional uniform gases with zero-range inter- 
action serve as a test ground for a wide class of one- 
dimensional systems. This was realized long ago by Lieb 
and Liniger, who used the Bethe-ansatz method @| to 
obtain the exact wave functions and spectrum of the one- 
component boson gas 0. They considered the following 
Hamiltonian for N identical bosons of mass m and chem- 
ical potential fi, 



H = - 



■ V — 



2m dxf 

i—i 



l<i<j<N 



5{x,-Xj)- tiN. (1) 



The interaction constant giD has the dimension of energy 
times length and is conventionally written in terms of a 
characteristic length any according to 



5lD 



2h^ 
maiD 



(2) 



The dimensionless coupling strength 7 involves the num- 
ber of particles per unit length, tiid = N/ L, where L is 
the length of the ID system, and is given by 



niDfliD 



(3) 



Thus when the density decreases, the coupling strength 
increases, and one approaches the Umit of iinpenetrable 
bosons, the so-called Tonks- Girardeau gas Q- In this 
limit the thermodynamic properties of the system are the 
same as those of a one-dimensional gas of non-interacting 
one-component fermions with the Fermi energy 
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Henceforth we set h = 2m — 1. These constants will be 
restored in Section IVll 

Similarly to the one-component case, one can consider 
a two-component gas of fermions (bosons). Such gases 
consist of particles with some internal quantum param- 
eter a, which can take two values. For brevity we refer 
to this parameter as "spin", and label the correspond- 
ing states by cr =tii • 111 the situation where the inter- 
action constant giu is tuned to be independent of the 
"spin" index, we take the first-quantized Hamiltonian of 
the system to be the same as for the one-component case, 
Eq. Q. The information about the system being two- 
component is encoded in the symmetry of the coordinate 
wave function. To introduce the "spin" indices in the 
Hamiltonian explicitly one can recast it in the second 
quantized form using the particle creation and annihila- 
tion operators t/^ti^) ^-nd '4'<t(x). Written in the second- 
quantized form the interaction term in Eq. Q is 



fflD 



for two-component fermions, and 



giD 
4 



dx : n{xY 



(5) 



(6) 



for two-component bosons. Here no- (a;) — ip^^{x)^pa{x) is 
the density of particles of type cr; the total density n{x) 
is n{x) = ni{x) + n.|(x), and the symbol :: stands for 
the normal ordering. For interactions of this form the 
total number of particles N and the total number N^^-^ 
of "spin" up (down) particles are good quantum numbers. 
Note, that the operator ^ is the most generic form of 
contact interaction for fermions, which is a consequence 
of the Pauli principle. In the bosonic case, however, the 
operator © represents a special situation where the three 
coupling constants are tuned to be equal. 

The density niD is controlled by the chemical potential 
fi, Eq. (QJ. In the infinite 7 limit the relation between 
niD and ^ does not depend on the statistics of particles 
and has the following simple form 



(7) 



For the one-component gas this result can be found in 
e.g. Refs. Q, for the two-component gas in Ref. Q. It 
follows from Eq. ^ that has a dimension of inverse 
length. Hereafter we measure distances in units of l/^/JI. 

What temperature scales characterize ID gases with 
the Hamiltonian (QJ in the limit of strong repulsion. 



7^1? For the one-component system there is only one 
temperature scale: the Fermi temperature Tp = E-p/kB, 
where fee is the Boltzmann constant. Low temperatures 
T 'ti Tp correspond to the degenerate quantum gas. 
The inclusion of "spin" gives rise to another character- 
istic temperature Tq <C Tp, defined by the bandwidth 
for the "spin" excitations. This bandwidth can be es- 
timated as follows. Consider the sector of the Hilbert 
space with a given number of particles N and a polariza- 
tion M = - Ni. Denote by E{M) the ground state 
energy of the Hamiltonian Q in this sector. For fermions 
the function E{M) has its absolute minimum at M = 
[l(l | in agreement with the Lieb-Mattis theorem • For 
bosons the function E{M) reaches its absolute minimum 
at \M\ = N giving rise to the demixing transition 
[T^ . The energy difference per particle between the po- 
larized (demixed) and unpolarized (mixed) phases 
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gives an estimate for the bandwidth of "spin" excitations: 
To ~ e/kB- In the large 7 limit e is easily found from the 
thermodynamic Bethe ansatz For fermions one finds 
that to the leading order in I/7 
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One can see that the characteristic temperature scale for 
"spin" excitations is given by Tq ~ Tf/t- The same esti- 
mate is valid for bosons. 

The two scales To and Tp are well separated when 7 
is large enough. In this situation one can distinguish 
between two different quantum (which means T ^ Tp) 
regimes: T < Tq and T > Tq. The first regime has been 
studied extensively and is conventionally described by 
the Luttinger Liquid (LL) theory (ijhthe results on the 
second one appeared very recentlv 11,4 IT^ for the 
limiting case Tq ^ T <C Tp and are qualitatively different 
from the LL theory predictions. This can be understood 
by noticing that for T > Tq "spin" degrees of freedom 
are strongly "disordered" thus violating the LL theory 
applicability conditions. At the same time, the density 
degrees of freedom are not affected by the temperature 
until it becomes of the order of Tp. 



III. ONE-PARTICLE DENSITY MATRIX: AN 
EXACT REPRESENTATION 

In this Section we investigate the one-particle density 
matrix of the two-component strongly repulsive (7 
+00) ID gases with the Hamiltonian Q for the two lim- 
iting cases T ^ To and To ^ T of the two quantum 
regimes discussed in Section |n| Our considerations will 
be based on the exact solution of the model Hamilto- 
nian Q in the infinite 7 limit. The two different regimes 
T <C Tq and Tq<^T are accessed by choosing the appro- 
priate order of limits T ^ and 7 ^ 00, as described in 



2 



[1^: To ensure Tq -^T one should take the limit 7 — > oo 
first and then set T = in the wave functions and spec- 
tra of the model, while for the LL regime T ^ Tq one 
should take the limit T ^ first. 

For two-component systems the one-particle density 
matrix is defined as 



(^|(x)^T(o)) + (V'|WV'i(o))' 



(10) 



Since the interaction terms Eqs. JHl and ^ are invariant 
with respect to "spin" rotation, one has (V'j(a;)V'T(0)) = 

(V'|(a;)V'i(0)). The density of particles is niD — 2p(0). 
Both physical intuition and crucial technical advance 
come from the observation that in the infinite 7 limit the 
exact Bethe-ansatz many-particle eigenfunctions of the 
Hamiltonian Q factorize into a coordinate part which is 
similar to the wave function of non-interacting spinless 
fermions and a spin part which corresponds to an eigen- 
state of the isotropic Heisenberg chain. Using this factor- 
ization Ogata and Shiba 0| derived an explicit formula 
for the density matrix lll()|l in the case of two-component 
fermions on a lattice (the Hubbard model). The two- 
component fermion gas being the Hmit of the Hubbard 
Hamiltonian for the vanishing fiUing factor, Ogata and 
Shiba's formula can be used. It also has a straightfor- 
ward extension to the bosonic case. Both bosonic and 
fermionic cases are discussed in the paragraph below. 

Using the results of Ref. we write the density ma- 
trix p{x) as 



(11) 



where the average is taken over the ground state of 
the system of non-interacting one-component fermions 
with creation and annihilation operators and x, re- 
spectively. The density of these fermions coincides with 
the density of physical particles tt-id- The operator 



dx'xHx')x{x') 



(12) 



counts the number of particles between points and x. 
Assuming lattice regularization with an infinitesimal lat- 
tice constant, Afx is a Hermitian operator with integer 
eigenvalues. To define the function uj{N) consider an 
infinite isotropic Heisenberg chain. Its Hilbert space is 
spanned by vectors \ . . . aicr2 ■ ■ . crj ■ . ■), where Cj =t, i is 
the spin at the j-th site. In this space one may intro- 
duce a string operator Qn, which acts as a cyclic shift 
operator on the string of spins of length N: 



r^Afl . . . (Jj(Jj + l ■ ■ ■ ^ j+N -2^ j+N -iC j+N ■ ■ ■) 
I • ■ -(yj+N-lCjCFj+l ■ ■ ■ <^j+N-2<^j+N ■ ■ 

The function uj{N) is defined as follows 



uj{N) 



N)H- 



(13) 



(14) 



Here (j) is a statistical angle: (j) = for bosons and 
= TT for fermions. In the case T <C To the average 
Q-H in Eq. ifnt is taken over the ground state of the 
isotropic Heisenberg chain, which is antiferromagnetic for 
the fermion gas and ferromagnetic for the boson gas. In 
the case T » Tq, the average is to be taken over all 
possible spin configurations. 

For the boson gas at exactly zero temperature and both 
fermion and boson gases at T 3> Tq the calculation of 
uj{N) is relatively simple. For bosons at zero tempera- 
ture the average in Eq. ifTiji is taken over the fully 
polarized ferromagnetic ground state of the spin chain 
which gives 



uj{N) 



(15) 



For bosons and fermions at T ^ Tq the average ( )-h in 
Eq. itTHl runs over all spin configurations. This gives 



u;{N) = e*('^-'^)^2 



-N 



(16) 



In both equations Eqs. II15|I and lllfill the function w(A^) 
is an exponential of TV : u!{N) ~ e'^^, for which the 
average ifTTll can be expressed in terms of the Fredholm 
determinant. More precisely, one obtains the following 
representation for the density matrix ljl()|l : 



p{x) 



P(0) 

Here the determinant 



det(/ + Vi + V2)ix) - det(/ + Vi)ix). (17) 
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V{kN,ki) 
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■■■ VikukN) 
■■■ V{kN,kN) 



(18) 



is the Fredholm determinant of a linear integral operator 
V with the kernel defined on [—1, 1] x [—1, 1]. In our case 
the kernels Vi and V2 are: 



Viik,p) 



sin |(fc — p) 
7r(fc — p) 
1 X 

V2{k,p) = -exp{-i-{k+p)} 



(19) 
(20) 



The values of the parameter a are given in the Ta- 
bleQI Recall that we measure distances in units of l/^/JI, 
Eq. 

It was mentioned above Eq. ((T5|l that the ground state 
of the two-component bosons is fully polarized in "spin" 
[1^ . Therefore, at T < Tq the bosons become effec- 
tively one-component (compare first and second lines 
in the Table ^ and we turn back to the original Lieb- 
Liniger model Q- The determinant representation Ijl7|l 
in this model was obtained by Lenard [l3|. For the two- 
component fermion gas at T Tq Eq. ifTTIl was obtained 
in Ref. p3|. Izergin and Pronko calculated both bosonic 
and fermionic density matrices in Ref. 
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a 


one-component boson gas, 


T = 


= 


-2 


two-component boson gas, 


T 


< 


To 


-2 


two-component boson gas, 


T 


> 


To 


-3/2 


two-component fermion gas. 


T 


< 


To 




two-component fermion gas. 


T 


> 


To 


-1/2 



Table L The values of the parameter a entering Eq. 11911 are 
given. For the two-component fermion gas at T <g To no 
exact results on p{x) are known. 



IV. ONE-PARTICLE DENSITY MATRIX: 
ASYMPTOTICS AND NUMERICS 

In this Section we analyze the short and long-distance 
asymptotics of the determinant representations of the 
density matrix (tTTHl . We also describe the numerical pro- 
cedure which deals with the intermediate length scales. 
The results of this Section will be used for the calculation 
of the momentum distribution function given in the next 
Section. 

Since p{x) = p{—x), we assume a; > 0. For short dis- 
tances (a; ^ 1) one can easily find from Eqs. II17|I and Ijl8|l 
that 



p{x) .T^ a x^ allx^ 

p(0) ^"^^y^^l8^120^ ^2700 



(21) 



For a = this is simply the expansion of sinx/x. No- 
tice that higher order terms that are not exhibited on the 
r.h.s. of Eq. (I21II . depend on a nonlinearly. Next, con- 
sider the long-distance expansion of Eq. IjlTII . Compared 
to the short-distance analysis, this is a more sophisticated 
task. For the boson gas at T <C Tq the analysis carried 
out in Refs. j23| shows the power-law decay of the density 
matrix: 



pjx) _ C 

p(o) Vi 



1 



1 



cos 2x 

4 



\Qx^ 



sin 2x 



(22) 

with relative corrections of the order of x ^. The con- 
stant C is given by C = 716^/22-1/3^-6 ^j^jj ^ ^ 

1.2824271 . . . being Glaisher's constant. For the two- 
component fermion gas at T » Tq the decay is expo- 
nential [H: 



p{x) = Cie-^'^^x'^^ sin(2a; - vlnx + C2). 



(23) 



Here Ci and are some constants which are given ex- 
plicitly in Ref. [22], while 



1/2 ln2 
Af---1, y = —- 

Z TT 



(24) 



The relative corrections to Eq. Ij23|l are of the order of 
x~^. The technique developed in the papers jlSvlGj and 
[T^ for the two-component fermion gas can be easily 
adopted to the case of the two-component boson gas at 



T ^ Tq. Like for fermions, the density matrix decays 
exponentially in the large x limit 



p{x) = const e"2''^a;'^= 
The anomalous exponent 

AB-i(^^-l) 



(25) 



(26) 



and v is given by Eq. ll24|l . The relative correction to 
Eq. lf25|l is of the order of x~^. 

Finally, we turn to intermediate values of x. To analyze 
the representation Ijl7|l . we use instead of Eq. Ijl8|l an 
alternative definition of the Fredholm determinant 23|. 
Let y be an LxL matrix with the entries Vab = V{ka, ki,). 
Let 



2a 
L - 1 



- 1 



a = 0,l,...,L-l. 



(27) 



Then the Fredholm determinant Eq. III8II can be repre- 
sented as follows: 



det(/ + V) = lim dctL / H 

L^oo \ L - 1 



V 



(28) 



The matrix / on the right hand side of the equation (2^ 
is the L X L identity matrix. We use the representation 
Ij28ll to calculate p{x) for x in the range < x < 20. For 
the figures shown in the next Section we used the number 
of divisions L — 800. 



V. LOW-TEMPERATURE RECONSTRUCTION 
OF THE MOMENTUM DISTRIBUTION 

In this Section we present our results on the momen- 
tum distribution function of quantum two-component 
gases in two temperature regimes T <C Tq and T ^ 
discussed in Section UTI 

The momentum distribution function n(fc) for a one- 
dimensional gas is given by 



n(fc) 



dxe 



-ikx 



(29) 



where p{x) is the one-particle density matrix ltT71|l . Recall 
that X is measured in units of l/y^ as defined below 
Eq. 0. 

First, consider the gas of two-component fermions. 
The momentum distribution for different temperature 
regimes is shown in Fig. The T <C Tq curve was ob- 
tained by the interpolation of the data given in Ref. 
for the Hubbard model at low fiUing. The T > To 
curve was calculated using the determinant representa- 
tion discussed in Sections IIIII and ITVl One can see a dra- 
matic change in the shape of the momentum distribution 
function as the temperature increases from T ^ Tq to 
T To . The singularity of the momentum distribution 
aX k = kp = 7rniD/2, predicted by the Luttinger theorem 



4 



1.0 

n{k) 

0.5 









■ 








1 





12 3 4 

k/kp 

Figure 1: The momentum distribution function n(fc) for two- 
component fermions. The thick solid curve is our result for 
T 2> To- The dashed curve corresponds to T <C To and is 
obtained by the interpolation of the data (filled circles) from 
Ref. The thin solid line shows the distribution function 

for non-interacting two-component fermions. The momentum 
k is normalized to the Fermi momentum ftp = 7miD/2. 

[23 |. disappears and the momentum distribution spreads 
out to larger k. Another very pronounced effect is that 
the large momentum tail of the momentum distribution 
gets strongly suppressed as the temperature increases. 
This is a counterintuitive result from the point of view 
of physics of weakly interacting systems, where the frac- 
tion of high energy particles grows monotonously with 
increasing temperature. In strongly correlated systems, 
however, bare particles have no overlap with elementary 
excitations of the system and the momentum distribution 
of bare particles is not directly related to the distribution 
of energy between the eigenmodes of the system. 

Next, consider the gas of two-component bosons. Us- 
ing the results of Sections IIIII and IIVI we calculate the 
momentum distribution function Eq. Ij29|l numerically for 
both T <C To and T ^ Tq. The results of the calculations 
are shown in Fig. [21 One can see that the momentum 
distribution in the ground state of bosons at T <C To 
is divergent at /c = 0. This behavior takes place at any 
coupling strength 7 (2R l2^ and is the manifestation of 
the quasi-long-range order in the system, that is, quasi- 
condensate. The situation is different at T To. Due 
to the exponentially decaying term in the asymptotic ex- 
pression Eq. Il25|l ■ the function n{k) is continuous with 
all its derivatives for all k. The suppression of the quasi- 
condensate component as the temperature increases from 
T ^ To to above Tq has a simple physical explanation: 
due to the excitation of the soft degrees of freedom the 
system loses coherence and cannot exhibit quasi-long- 
range order. For the large momentum tail of n{k) one 
gets from Eq. Ij21ll : 

nik) ^, k^oo. (30) 

k^ 

Comparing the values of a given in Table U one can see 
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Figure 2: The momentum distribution function n{k) for two- 
component bosons. The solid curve is the result for T ^ To, 
the dashed curve is for T ^ To . For convenience we normalize 
the momentum k to the Fermi momentum fcp — nnio /2 of a 
fermion system with the same density. 

that the fraction of particles with large momentum de- 
creases as the temperature increases from T ^ To to 
above Tq similarly to the fermionic case. 

VI. DISCUSSION 

One-dimensional quantum gases were recently created 
and studied in an optical lattice forming an array of well 
separated one-dimensional tubes Q. The particle mo- 
mentum distribution in such systems can be observed in 
experiments involving ballistic expansion |^ . Large 
values of the coupling constant, 7 ~ 200, were reported 
in Ref. Such values of 7 would give a two orders 
of magnitude separation between the two temperature 
scales To and Tp, which is sufficient for the momentum 
reconstruction to be observable. There are, however, cer- 
tain geometric limitations on the validity of the model, 
Eq. 1^, considered in this paper, which we discuss be- 
low. For the convenience of the reader we restore the 
constants h and m in the equations. 

Consider an elongated harmonic trap of transverse fre- 
quency uj± and axial frequency ujz = Awj,, A ^ 1, 
containing a quantum gas of N particles. The char- 
acteristic length aiD in Eq. Q is related to the three- 
dimensi onal sca ttering length a and the transverse width 
a± — ^JKfmu^ by (23 

aiD-^fl-Ci — V (31) 

where Ci = -C(l/2)/\/2 = 1.0326 . . . is a numerical con- 
stant, and C^{x) is the Riemann zeta function. Since 
the local chemical potential /^(z) is given by ^i{z) = 
7r^ft^n^j-,(z)/2m, the one-dimensional density is found 
within the Thomas-Fermi approximation to be nij:,{z) = 
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{2N — z"^ / alY^'^ /iraz, where = yjh/mujz- The value of 
7, Eq. may then be expressed in terms of the particle 
number TV, the scattering length a and the trap frequen- 
cies [jj_L and ujz. At the trap center (z = 0) one obtains 
the simple expression 

From Eq. ll32|l it is clear that to achieve large values of 
7(0) at fixed TV and A the system must be tuned to the 
resonance — Cia = 0. However, very close to the reso- 
nance the scattering amplitude will no longer correspond 
to a delta-function potential in real space due to the pres- 
ence of higher-order terms in the relative momentum. 
This may be demonstrated from the expansion of the 
scattering amplitude / as given in [2^ . 



l + ikaiiy-iV2C{i/2){kai_Y/%' 

From this expression one notes that the terms of 
third order in k are negligible only when aiD/a_L ^ 
y2C(3/2)(fca±)2/8. If we identify k with the Fermi 
wave number fcp which has its maximum value fcp = 
7miD(0) = ^/2N /az in the center {z ~ 0) we obtain 
the condition 

a,^yya^^^^^NX (34) 

which in combination with Eq. (13 111 implies that (aj^/a — 
Ci)/N\ y> 1, in order for our model to be appHcable. 



Combining this result with Eq. (I32II we arrive at the up- 
per bound for the interaction constant 7 for the trap with 
given A and N: 

3 

1 « 7(0) « 7tV2 (^-^^ ' . (35) 

For traps with A ^ 300 containing iV = 30 particles 
this estimate gives 7 <C 140. 



VII. CONCLUSIONS 

In conclusion, the momentum distribution of two com- 
ponent boson and fermion gases was considered in the 
limit of strong interatomic repulsion. It was shown that 
due to the strongly correlated nature of the system the 
presence of the internal degrees of freedom of the atoms 
results in the emergence of a new low temperature deco- 
herent state absent in one-component gases. The onset 
of this state with increasing temperature is marked by a 
pronounced change of the momentum distribution of the 
atoms. 

Acknowledgement 

The authors would Hke to thank J. DaHbard for helpful 
discussions. M.B. Zvonarev's work was supported by the 
Danish Technical Research Council via the Framework 
Programme on Superconductivity. 



[1] A. Gorlitz et al, Phys. Rev. Lett. 87, 130402 (2001); 
F. Schreck, ibid. 87, 080403 (2001); S. Dettmer, ibid. 87, 
160406 (2001); 

[2] H. Moritz, T. Stoferle, M. Kohl, and T. Esslinger, Phys. 
Rev. Lett. 91, 250402 (2003); B. Laburthe Tolra et ai, 
Phys. Rev. Lett. 92, 190401 (2004). 

[3] S. Richard, Phys. Rev. Lett. 91, 010405 (2003). 

[4] B. Paredes et al, Nature 429, 277 (2004). 

[5] M. Kohl et al, e-print: cond-mat/0406397 

[6] V.E. Korepin, N.M. Bogoliubov and A.G. Izergin, Quan- 
tum Inverse Scattering Method and Correlation Func- 
tions (Cambridge University Press, Cambridge, England, 
1993). 

[7] E.H. Lieb and W. Liniger, Phys. Rev. 130, 1605 (1963); 

E.H. Lieb. ibid. 130, 1616 (1963). 
[8] L. Tonks, Phys. Rev. 50, 955 (1936); M. Girardeau, .1. 

Math. Phys. (N.Y.) 1, 516 (1960). 
[9] A.G. Izergin and A.G. Pronko, Nucl. Phys. B 520, 594 

(1998). 

[10] C.N. Yang, Phys. Rev. Lett. 19, 1312 (1967). 
[11] E. Lieb and D. Mattis, Phys. Rev. 125, 164 (1962). 
[12] You-Quan Li et al, Europhys. Lett. 61, 368 (2003). 
[13] M.A. Cazalilla and A.F. Ho, Phys. Rev. Lett. 91, 150403 
(2003). 

[14] A.O. Gogolin, A. A. Nersesyan, and A.M. Tsvelik, 



Bosonization and Strongly Correlated Systems (Cam- 
bridge University Press, Cambridge, England, 1998). 

[15] Vadim V. Cheianov and M.B. Zvonarev, Phys. Rev. Lett. 
92, 176401 (2004). 

[16] Vadim V. Cheianov and M.B. Zvonarev, J. Phys. A: 
Math. Gen. 37, 2261 (2004). 

[17] G.A. Fiete and L. Balents, e-print: Lcond-mat/0403744j 

[18] M. Ogata and H. Shiba, Phys. Rev. B 41 2326 (1990). 

[19] A. Lenard, J. Math. Phys. (N.Y.) 5, 930 (1964). 

[20] A. Berkovich and J.H. Lowenstein, Nucl. Phys. B 285, 
70 (1987). 

[21] K.B. Efetov and A.I. Larkin, Sov. Phys. JETP 42, 390 

(1976); H. G. Vaidya and C. A. Tracy, Phys. Rev. Lett. 

42, 3 (1979); M. Jimbo, T. Miwa, Y. Mori, and M. Sato, 

Physica (Amsterdam) ID, 80 (1980). 
[22] A. Berkovich, J. Phys. A 24, 1543 (1991). 
[23] V.I. Smirnov, A course of higher mathematics, Vol IV, p. 

24 (Pergamon, Oxford, 1964). 
[24] J.M. Luttinger, Phys. Rev. 119, 1153 (1960); K.B. Bla- 

goev and K.S. Bedell, Phys. Rev. Lett. 79, 1106 (1997). 
[25] V.N. Popov, JETP Lett, 31, 526 (1980). 
[26] F.D.M. Haldane, Phys. Rev. Lett. 47, 1840 (1981). 
[27] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 



6 



